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Abstract 

We investigate numerically the three-dimensional 0(2) model on 8 3 -160 3 lattices as 
a function of the magnetic field H. In the low-temperature phase we verify the in- 
dependence of the magnetization M induced by Goldstone modes and determine M 
for V — > oo on the coexistence line both by extrapolation and by chiral perturbation 
theory. This enables us to calculate the corresponding critical amplitude. At T c the 
critical scaling behaviour of the magnetization as a function of H is used to determine 
another critical amplitude. In both cases we find negative corrections-to-scaling. 
Our low-temperature results are well described by the perturbative form of the 
model's magnetic equation of state, with coefficients determined nonperturbatively 
from our data. The 0(2) scaling function for the magnetization is found to have a 
smaller slope than the one for the 0(4) model. 
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1 Introduction 



0(N) models are of general relevance to condensed matter physics and to quantum 
field theory, because many physical systems exhibit a second-order phase transition 
with universal properties, which may be derived analytically from O(N) symmetric 
effective field theories |TJ. A further aspect of these theories is remarkable. Due 
to the existence of massless Goldstone modes in O(N) models with N > 1 and 
dimension 2 < d < 4 [§, [3| singularities are expected on the whole coexistence line 
T < T c , H = 0, in addition to the known critical behaviour at T c . For N = 2 and 
d = 3 there exists a rigorous proof f4|. Numerically these theoretical predictions 
have been recently confirmed by simulations of the 3d 0(4) model ||. In the same 
article a parametrization of the equation of state including the Goldstone effect has 
been worked out and compared to perturbative predictions [H, [/J. Here we want 
to extend these numerical studies to the case of the 0(2) model. For the latter 
approximate representations of the equation of state have been derived, starting 
from high-temperature expansions ||. 

Our interest in the 0(2) model is moreover motivated by its relation to the 
staggered formulation of quantum chromodynamics (QCD) on the lattice. At finite 
temperature QCD undergoes a chiral phase transition. For two degenerate light- 
quark flavors this transition is supposed to be of second order in the continuum 
limit and to belong to the same universality class as the 3d 0(4) model |Q-[TT|. 



QCD lattice data have therefore been compared to the 0(4) scaling function, as 
determined numerically in fl2 |. For staggered fermions the comparison to 0(4) is 



still not conclusive fI5|-||15|) but results for Wilson fermions |]T6| seem to agree quite 
well with the predictions, though for the Wilson action the chiral symmetry is only 
restored in the continuum limit. In the staggered formulation a part of the chiral 
symmetry is remaining even for finite lattice spacing, and that is 0(2). For the test 
of lattice QCD with staggered fermions it is therefore also important to know the 
corresponding universal 0(2) scaling function or equation of state. 

The 0(2)-invariant nonlinear a- model (or XY model), which we want to inves- 
tigate in the following is defined by 

pn = -j J2 s *- s ; - H E s *> (!) 

<i,3> i 

where % and j are nearest-neighbour sites on a d— dimensional hypercubic lattice, 
and Sj is an 2-component unit vector at site i. It is convenient to decompose the 
spin vector S« into a longitudinal (parallel to the magnetic field H) and a transverse 
component 

s, = sfn + si . (2) 

The order parameter of the system, the magnetization M, is then the expectation 
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value of the lattice average 5" of the longitudinal spin component 

M =<^E^I > = < 5 "> • (3) 

i 

On finite lattices and H = system flips occur and lead to < 5" >= 0. Therefore 
one usually resorts to approximate order parameter definitions, as e.g. < \S"\ > 



- see the discussion in Tapalov and Blote fll7| . Here this is unnecessary since we 
always work at finite H . 

There are two types of susceptibilities. The longitudinal susceptibility is defined 
as usual by the derivative of the magnetization, whereas the transverse susceptibil- 
ity corresponds to the fluctuation of the lattice average S -1 of the transverse spin 
component 

XL = H = V« S" 2 > — M 2 ) , (4) 
XT = V<S ±2 > . (5) 



Both susceptibilities are predicted to diverge on the coexistence line T < T c , H = 0. 
The Goldstone singularities lead for all temperatures below T c for small H to strong 
finite-size effects, which have been studied in the context of chiral perturbation 
theory for the three-dimensional 0(2) model already in JTE[] . An explicit check on 
the if-dependence of the magnetization for V — > oo close to H = is however 
lacking. In comparison to the 0(4) model, we expect the Goldstone effects to be 
weaker, because we have only one transverse spin component in 0(2), not three like 
in 0(4). Another difference to 0(4) is that sizeable corrections to scaling appear 
in the XY model ||19|| . The determination of the universal equation of state of 
the three-dimensional 0(2) model, therefore requires a careful consideration of the 
corrections. We could have avoided the problems with these corrections by working, 



as it was done in Ref. |19j at A = 2.1, instead of A = oo. For a better comparison 
to the 0(4) case ||, we chose however to use the same model definition. Also there 
is then no need for an additional modulus update. 

The plan of the paper is as follows. In the next section we briefly discuss the 
perturbative predictions for the magnetization and the susceptibilities at low tem- 
peratures, the form and the analytic results for the magnetic equation of state. A 
more general review can be found in ||. Our numerical results are presented in 
Section [5], the equation of state is determined and parametrized in Section |j. We 
close with a summary and our conclusions in Section [5|. 



2 Perturbative Predictions and Equation of State 

For T < T c the system is in a broken phase, i.e. the magnetization M(T, H) attains a 
finite value M(T, 0) at H = 0. We explicitly assume here H > 0. As a consequence 
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the transverse susceptibility, which is directly related to the fluctuation of the Gold- 
stone modes, diverges as H~ l when H —>■ for all T < T c . This is immediately clear 
from the identity |J 

M(T, H) 

Xt = jj , (6) 

which is valid for all values of T and H. It is non-trivial that also the longitudinal 
susceptibility is diverging on the coexistence line for 2 < d < 4. The leading term in 
the perturbative expansion for 2 < d < 4 is H d l 2 ~ 2 |2|, [5J. The predicted divergence 
in d = 3 is thus 

X l(T<T c ,H) ~ H- 1 ' 2 . (7) 

This is equivalent to an if 1,/2 -behaviour of the magnetization near the coexistence 
curve 

M(T < T C ,H) = M(T,0) + c H 1/2 . (8) 

An interesting question is whether this form of the magnetization is compatible with 
the general Widom-Griffiths form of the equation of state describing the critical 
behaviour in the vicinity of T c . It is given by 

y = f(x) , (9) 

where 

y = h/M s , x = t/M 1/p . (10) 

The variables t and h are the normalized reduced temperature t = (T — T c )/To and 
magnetic field h = H/H . We take the usual normalization conditions 

/(0) = 1, /(-1) = 0. (11) 

The critical exponents 5 and (3 appearing in Eqs. ^ and |1(] specify all the other 
critical exponents 

dv = P(l + 6), 7 = /3(5-l), v c = v/(35. (12) 

Possible dependencies on irrelevant scaling fields and exponents are however not 
taken into account in Eq. ^, the function f(x) is universal. Another way to express 
the dependence of the magnetization on t and h is 

M = hV'fait/h 1 '*) , (13) 

where fo is a scaling function. This type of scaling equation is used for comparison 
to QCD lattice data. The scaling forms in Eqs. (0) and fli"3]) are clearly equivalent. 



The equation of state (0) has been derived by Brezin et al. to order e 2 in the 
e-expansion, where e = 4 — d. The resulting approximation has been considered by 
Wallace and Zia [0] in the limit x — > — 1, i.e. at T < T c and close to the coexistence 
curve. In this limit the function is inverted to give x + 1 as a double expansion in 
powers of y and y d / 2 ~ l 

x + 1 = c lV + d 2 y d/2 - 1 + d lV 2 + d 2 y d / 2 + d 3 y d - 2 + ... . (14) 
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The coefficients ci, £2 and d% are then obtained from the general expression of ||. 
The above form is motivated by the if -dependence in the e-expansion of xl at low 
temperatures f7|. In Section £| we propose a fit of our Monte Carlo data to the 
perturbative form of the equation of state, using (OX 



As for the large-x limit (corresponding to T > T c and small if), the expected 



behaviour is given by Griffiths's analyticity condition [20 



f(x) = ^2a n x^-W. (15) 

n=l 



3 Numerical Results 



Our simulations are done on three-dimensional lattices with periodic boundary con- 



ditions and linear extensions L = 8 — 160 using the cluster algorithm of Ref. |2"I 
The value J c = 0.454165(4), obtained in simulations of the zero-field model [22], 
is used in the following. We have simulated at increasingly larger values of L at 
fixed values of J = 1/T (i.e. at fixed temperature T) and if in order to eliminate 
finite-size effects. 

As an example we show in Fig. [1] part of our data for the magnetization for 
low temperatures, at J = 0.5, 0.47, and at J c , plotted versus if 1 / 2 . The picture 
is rather similar to the one obtained in 0(4) ||: strong finite-size effects appear 
for small if and persist as one moves away from T c , the results from the largest 
lattices are at first sight linear in if 1//2 , as predicted by Eq. |8], below but not at T c . 
A closer look at the curves reveals however subtle differences compared to the 0(4) 
case, which are confirmed by the attempt to determine the value of M(T, H = 0) 
by a linear fit in if 1 / 2 . Very close to if = the fixed temperature curves become 
slighty flatter, leading to a higher value for M(T, 0) than expected from the data 
at larger if values. This behaviour is more pronounced close to T c than at lower 
temperatures. The determination of the magnetization in the thermodynamic limit 
(V — > 00) on the coexistence line requires therefore more precise data and smaller 
values of H than for the 0(4) model. In order to extrapolate the data to H —>■ 
and V — > 00 we apply two different strategies. The first is to extend the linear form 
in if 1 / 2 , Eq. [|, to a quadratic one @] 

M(T < T c , H) = M(T, 0) + ci if 1/2 + c 2 ff , (16) 

and to fit the data from the largest lattices, which we assume to represent data 
on an infinite volume lattice, to this form. The second way to find M (T, 0) is just 
opposite to the first. Here we exploit the L or volume dependence at fixed J and 
fixed small if to determine via chiral perturbation theory (CPT) the magnetization 
S of the continuum theory for V — > 00, ff = 0, which is related to M (T, 0) by 
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Figure 1: The magnetization as a function of H 1 ^ 2 in the low-temperature region 
for fixed J = 0.5, 0.47 and J c , starting with the highest curve for various L. 



m(t ' 0) = 7j- 



(17) 



According to the e-expansion of CPT the quantity £ and the Goldstone boson decay 
constant F are the only parameters which determine the finite-size corrections to 
the partition function to order 1/L 2 [E3]. We summarize the relevant formulae for 



the three-dimensional XY model from Ref. |I8 |. In our notation the magnetization 
is 



M 



w[pi?7 + 2p 2 a 2 



where u = p^HV/VJ, a = 1/(F 2 L) and 

Pi = l+ l -{3ia + ^{(3l-2{3 2 )a 2 



P2 = -^2 



(19) 



For a symmetric three-dimensional box j3i = 0.225785 and (5 2 = 0.010608. The 
parameter rj is given in terms of modified Bessel functions as 



n 



h{u) 
ul 2 (u) 



(20) 



By construction the e-expansion is only applicable in a range where m n L<l which 
translates into the condition 



E f F s 2 



(21) 
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In Fig. we show the data for the magnetization as a function of H 1 ^ 2 for six fixed 
values of J in the low-temperature phase. We observe a remarkable coincidence of 
the fits according to Eq. |16| with the CPT results at H = 0. Details of the fits are 
presented in Table |l|. The y 2 per degree of freedom for the first fit type is in the 
range 0.01-0.5 . The corresponding numbers for the CPT fits are of the order of 1. 
Like for the 0(4) model we see in Fig. |] that the region where M is linear in H 1 ^ 2 
shrinks when T approaches T c , the value of c 2 in Table [l] increases correspondingly. 
At J = 0.460 the fit parabola coincides only on a small piece with the largest volume 
data (here we simulated even on L = 160 lattices), the errors on c\ and c 2 become 
large, though the value at H = agrees again nicely with the CPT result. We have 
further convinced ourselves that M(T, 0) can also be obtained from calculations at 
H = with the modulus definition. The approach to the thermodynamic limit is 
here from above. 
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Figure 2: The magnetization vs. H 1 / 2 in the low-temperature region for fixed J = 
0.47, 0.465, 0.462 (left plot) and J = 0.6, 0.55, 0.5 (right plot) and different L 
(notation as in Fig.|I]). The filled circles are the results from CPT, the dashed lines 
the fits according to Eq. [16|. The solid lines connecting the data are from reweighting. 



In the neighbourhood of the critical temperature the results for M(T, 0) should show 
the usual critical behaviour. Contrary to the 0(4) case we expect here a sizeable 
correction to the leading scaling behaviour In order to determine the critical 
amplitude B we therefore make the following ansatz 

M(T < T c , H = 0) = B{T C - Tf[l + b x {T c - T) wv + b 2 (T c - T)} . (22) 



As an input we take the critical exponents from Ref. ]T9 



(3 = 0.3490(6) , v = 0.6723(11) , u = 0.79(2) . (23) 
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Table 1: Parameters of the fit of M(T, H) to Eq. [16] in the if -range of column 5 and 
the results of the CPT fit for H = 0.0001 in the L-range of column 7. 



A fit of all points of Table |], apart from the one at J = 0.46, using the form ([22]) 
leads to the result B = 0.945(5), b x = -0.053(23) and b 2 = -0.098(23). In Fig. get 
we show this fit and also the leading term separately. 
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Figure 3: (a) The magnetization on the coexistence line versus T — T c , the points 
are from Table [l], the solid line is the fit (0), the dashed line its leading part, 
(b) The magnetization at T c as a function of H, the line is the fit (|2"4|), the dashed 
line the leading part. 



As the critical point is reached the if-dependence of the magnetization changes to 
satisfy critical scaling. We therefore fit the data from the largest lattice sizes at T c 
to the form 

M{T C ,H) = d c H 1/s [l + d 1 c H U}Uc ] . (24) 
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A further term proportional to H is unnecessary, because the corrections to scaling 
are much smaller here, than on the coexistence line. The largest L data can be fitted 
very well with the ansatz ( p4|) and the critical exponents as input, as can be seen in 
Fig. |b. We find d c = 0.978(2) and d\ = -0.075(5). 

We have performed in addition some simulations in the high temperature phase of 
the model. Here the finite size effects are rather small. In Fig. f| we show the results. 
One observes that with increasing T (decreasing J) the region of linear dependence 
on H of the magnetization becomes larger. 
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Figure 4: The magnetization as a function of H in the high-temperature region for 
fixed J = 0.44, 0.4 and 0.3, starting with the highest curve for various L. 



4 Determination of the Equation of State 

In the last section we derived from our Monte Carlo data the magnetization in the 
thermodynamic limit. At temperatures below T c this was achieved for the larger 
H values by simply performing simulations on lattices with increasing size until all 
volume dependence was gone. For the small H region down to H = we took then 
advantage of the Goldstone effect to extrapolate our data and confirmed this by the 
CPT results for E. In the high-temperature phase we reached the thermodynamic 
limit already on lattices with L < 72 for all H values. 

The two critical amplitudes B and d c can now be utilized to normalize temperature 
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and magnetic field according to Eq. ([TT|). We find from our fits ( p2|) and (p^) 

T = B~ 1/p = 1.18(2) , and H = d' s = 1.11(1) . (25) 

As mentioned already, the equation of state does not account for possible corrections 
to scaling. A more general form of Eq. (|13|) is 

Mhr 1/S = V(ttr 1/0S ,h UVc ) . (26) 

Expanding the function \1/ in h UJUc leads to 

Mh~ l/S = f G {th- 1/p5 ) + h UJ ^f { G ) {th- 1/p5 ) + --- . (27) 

In order to obtain the scaling function fa we therefore perform quadratic fits to 
our data in h UUc at fixed values of th~ 1 ^ 5 in the low-temperature region, where 
the corrections are strong. Fortunately, the corrections are very small in the high- 
temperature region and the data scale directly. Altogether data with H < 0.0075 
and 0.43 < J < 0.55 were used for this purpose. In Fig. |5] we have plotted the 
uncorrected and the final results for the scaling function fa- 



2 



l 





ii 




II 










1 . 5 


i 


V 

\ 




1 . 


i 








3-2-1 


I I 



-20 -15 -10 -5 



Figure 5: The scaling function fg, Eq. (|i3|), (solid line). Also shown are the results 
for M/h 1 / 5 at fixed values of J = 0.55, 0.50, 0.47, 0.465 and 0.462 (dashed lines) 
which were used for the extrapolation, starting with the lowest curve. The circles 
are single data points. The inset shows the low-temperature region close to T c , with 
data for J < 0.46 (circles), 0.462 (triangles), 0.465 (crosses) and 0.47 (squares). 

Like in the case of 0(4) [|5| we want to describe the low-temperature part of the 
equation of state by its perturbative form as discussed in Section §, but with non- 
perturbative coefficients. We then want, again as in Ref. ||, to interpolate this 
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result with a fit to the large-x form (|T5|). The variables x and y are simply related 
to the scaling function Jq and its argument by 

y = fa 5 , x = (t/h^ s ) r G llP • (28) 



We first perform a fit to x in the interval [—1, 1.4] (which corresponds to y G [0, 3]) 
with the three leading terms in (IT3) 

xi(y) + 1 = (ci + d 3 ) y + c 22 / 1/2 + d 2 ?/ 3/2 . (29) 

In the fit we require di = 1 — (ci + c?3 + c 2 ) to fix the normalization y(0) = 1. We 
obtain 

c x + d 3 = 0.352(30) , c 2 = 0.592(10) . (30) 

The fit describes the corrected scaling function at T < T c and and the direct data 
in the high temperature phase up to x ~ 1.7. This confirms, as in 0(4), that the 
expression (14) is valid also away from x ~ — 1. Our coefficients can be compared 
to those calculated perturbatively for iV = 2 in Ref. [0] 

cj + d 3 = 0.818 , c 2 = 0.229 . (31) 

The result ( |30D is much closer to the coefficients found in 0(4) |5] than to the e- 
expansion results (|31]), though the Goldstone effect is decreasing somewhat with 
decreasing N: c 2 and the x region where the ansatz (|29|) is valid are a little smaller. 




Figure 6: (a) The function y = f(x) from fits at small x (solid line) and at large x 
(dashed line), (b) The interpolation (|34D for f(x) (solid line), and the parametriza- 
tions (A) (short dashs) and (B) (long dashs), n = 1 of f(x) from Ref. 0. In both 
parts we show high-temperature data ( filled circles). 
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For large x we have done a 2-parameter fit of the behaviour (]I5D, in the corresponding 
form for x in terms of y 

x 2 (y) = ay lh + by^ 2 ^ . (32) 
Considering data points with y in the interval [4,23000] we find 

a = 1.2595(30) , b = -1.163(20) . (33) 



Expression ( |32"D describes the data for x > 1.5. The small- and large-x curves cover 
the whole range of values of x remarkably well. The two curves overlap around 
y fa 3.5. This is shown in Fig. |6|a. We therefore interpolate them smoothly with the 
ansatz g 

x{y) = x x {y) 6 V ° + x 2 {y) 6 V 6 , (34) 

Do + y b vo + y 

where y$ = 3.5. In Fig. ^|b we compare our interpolation (|34]) of the equation of 
state to the high-temperature data and to two parametric forms for f(x) obtained 



from the high-temperature expansion in Ref. ||. According to the authors |24 
the difference between these two curves gives an idea of the uncertainty of their 
approach. For large x this difference can be traced back to the uncertainty in the 
universal amplitude ratio R x (Eq. (49) of Ref. ||), which was determined in || to 
R x = 1.4(1), that is with a 7% error. In fact, R x is related to the quantity a of our 
large- x parametrization by 

R x = a 7 = 1.356(4) , (35) 

which is in agreement with the result from Ref. || . Their low temperature parameter 
Cf may be calculated from our value of c 2 and amounts to c/ = 2.85(7). Finally 
we show in Fig. |7| the scaling function fa obtained parametrically from x(y) in d34|) 
and the corresponding scaling function for 0(4). As can be seen from Fig. 0, the 
dependence on t/h 1 ^ 13 is similar, but the 0(2) curve is flatter than the one for 0(4). 



5 Summary and Conclusions 

We have simulated the three-dimensional 0(2) model on cubic lattices as a function 
of the magnetic field H and the temperature T. From the behaviour of the measured 
magnetization M below the critical temperature T c and close to the coexistence line 
(H = 0), we could clearly verify the Goldstone-mode effects. This was done in 
two independent ways which led to the same values M(T < T c , H —> 0) for the 
magnetization in the thermodynamic limit V — > oo : on one hand we were able to 
profit by the observed Goldstone behaviour to extrapolate our data to H — > and 
on the other hand this result was confirmed from the finite-size dependence induced 
by the Goldstone modes, using chiral perturbation theory. 
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Figure 7: The scaling function fo = M/h 1 ^ 5 for the 0(2) model (solid line) and the 
0(4) model (dashed line). 

The M values on the coexistence line were subsequently used to calculate the 
critical amplitude B of the magnetization. Not unexpectedly, we found here strong 
corrections to scaling. On the critical line T = T c we computed then the second 
critical amplitude d c of the magnetization. Here, the corrections to scaling are less 
pronounced ; in the high-temperature phase they are unimportant. 

The problem with corrrections to scaling appeared again in the determination 
of the equation of state in the low-temperature phase. By generalizing the scaling 
equation to include possible corrections, we were able to derive the universal form of 
the equation of state for the 0(2) model from our data. This form was then analyzed 
in a similar manner as it was done for the 0(4) model in Ref. ||. In particular, we 
found again an efficient parametrization of the equation of state in the low x region 
with the perturbative form (^), which is based on expansion ([TJ|) by Wallace and 
Zia 0. The effect due to the Goldstone modes is nearly as large as in the 0(4) case, 
as can be seen from the coefficient £2 and much larger than the prediction from the 
e-expansion. 

Like for 0(4) we found a very good large-x fit to the high-temperature data. The 
coefficient a obtained from the fit implies a value of R x = 1.356(4), in agreement 
with parametrizations of the equation of state by Campostrini et al. ||. A further, 
indirect check of the equation of state is the computation of the universal ratio 
A + /A^. This requires however an integration of the magnetization with respect to 
h and derivatives with respect to t. We shall consider this in the near future. 
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Upon interpolation with the low-x curve, a complete description for the equa- 
tion of state is obtained, which can be plotted parametrically also for the scaling 
function. In Ref. p[ the lattice QCD data of the MILC Collaboration for iV r = 4 
were compared to the 0(4) scaling function. The test failed because the data were 
indicating a steeper scaling function. Since the 0(2) scaling function is flatter than 
the one for 0(4), the situation will be worse there. A way out may be the compar- 
ison to finite-size-scaling functions, since lattice QCD is presumably far from the 
thermodynamic limit. We are investigating this currently. 
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